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Abstract 


Lemon verbena leave is a flavoring food additive as well as a good source of valuable compounds such as 
essential oils, flavonoids and phenolic acids. However, similar to many other aromatic plants, lemon verbena 
leave is perishable due to its high moisture content. The aim of this work was to study the effect of air 
temperature (45, 55, and 65°C) on the quality attributes of lemon verbena leaves during hot-air drying (HAD). 
The drying kinetics were also modeled. The results showed that higher drying temperature led to a significant 
decrease (p<0.05) in the rehydration ratio due to a change in the structural features of the dried leaves. The 
essential oil content of dried samples was also significantly different (p<0.05) from that of the fresh leaves due 
to high loss of volatile components and ranged from 0.42 to 0.85. Moreover, a significant increase in the value of 
effective moisture diffusivity (Dem) and color change was observed when the samples were dried at 65°C 
compared to 45°C. The value of Derr varied from 1.140x10°!° to 2.280x10° m7/s and the activation energy was 
found to be 31.04 kJ/mol. The greatest R*? (©0.999) and the lowest RMSE and SSE were obtained for the 
Naghavi et al. model (proposed in this research) 
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Introduction 


Lemon verbena (Lippia citriodora) is a type 
of herb which is widely raised in western 
South America. It is also cultivated in Iran and 
mainly consumed as a spice and a medicinal 
plant (Funes ef al., 2009). There is an 
increasing interest in using lemon verbena 
leave in the food industry, because it is 
generally considered as a flavoring food 
additive. The leaves of lemon verbena have 
compounds such as essential oils, flavonoids 
and phenolic acids, which possess antioxidant 
activity (Pereira et al., 2007). They are mainly 
used to make herbal teas and refreshing 
sorbets as well as creating a lemon flavor in a 
number of food products such as fish and 
poultry dishes, jams, salad dressings, 
puddings, and beverages (Funes et al., 2009). 
Moreover, the leaves have digestive, sedative, 
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antispasmodic, stomachic, and antipyretic 
properties (Pereira et al., 2007; Funes et al., 
2009). However, similar to many other 
aromatic plants and herbs, lemon verbena 
leaves are perishable to microbial growth, 
mainly due to their high moisture content 
(around 84-85% wet basis). 

Drying is used to extend the shelf life of 
fruits, vegetables and aromatic plants as well 
as to reduce or suppress their enzymatic and 
microbial activities (Doymaz 2009; Doymaz 
2012). Aromatic plants are dried in order to 
extract their valuable compounds by solvents. 
Among the drying methods, hot-air drying 
(HAD) is still the most popular method, which 
is being employed to decrease the moisture 
content of foods and plants. Although HAD is 
time and energy consuming (Erbay and Icier 
2010), it has gained considerable attention by 
researchers due to its low capital cost 
compared to other drying techniques such as 
freeze-drying and infrared-drying. For this 
reason, it is still extensively employed by 
many researchers for long-term preservation of 
foods and herbs (Erbay and Icier 2010; 
Doymaz 2012; Lemus-Mondaca et al., 2015; 
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Oberoi and Sogi 2015; Aral and Bese 2016; 
Nozad et al., 2016; Roshanak ef al., 2016; 
Salarikia et al., 2016). 

Effective moisture diffusivity (Der) and 
activation energy (Ea) are two important 
physical properties of dehydrated foods which 
represent the rate of moisture loss during HAD 
and the level of energy needed to initiate a 
chemical reaction and to activate moisture 
diffusion, respectively (Lemus-Mondaca et al., 
2015). Derr can highly affect the drying 
kinetics and consequently the quality of 
dehydrated foodstuffs. Also, mathematical 
modeling of mass transfer during HAD 
requires the values of Derr. On the other hand, 
rehydration ratio (RR) is a key physical 
characteristic of dried foods which can reflect 
the degree of textural damage (such as 
shrinkage and tissue collapse) to foodstuffs 
during HAD (Doymaz et al., 2015). 

Numerous research papers can be found on 
the determination of Der, Ea, and RR during 
HAD of various aromatic plants and herbs 
(Doymaz 2012; Tasirin et al., 2014; 
Lemus-Mondaca et al., 2015; Nozad et al., 
2016; Salarikia et al., 2016). However, to the 
best of our knowledge, no study has been 
conducted on the determination of Des, Ea, 
RR, and color change for lemon verbena 
leaves under HAD conditions. The purpose of 
this research was to investigate the influence 
of hot-air temperature on the drying kinetics, 
color change, RR, and essential oil content of 
lemon verbena leaves under HAD and to 
calculate Derr and Ea, as well as empirical 
modeling of the dimensionless moisture ratio 
as a function of drying time. 


Materials and methods 
Materials 


Fresh lemon verbena leaves were collected 
every morning from a farm located in Tabriz 
(Iran), and immediately transferred to the 
laboratory. The leaves were sorted visually 
based on size, shape, color, and freshness and 
stored under refrigerated conditions (at 5°C) 
(Lemus-Mondaca et al., 2015) until use. The 
initial moisture content of the leaves was 


measured using the AOAC method (AOAC 
1984) and found to be equal to 84.72% (wet 
basis). 


Hot-air drying 


First, the leaves were removed from the 
refrigerator and arranged uniformly as a thin 
layer in a stainless steel basket. Then, they 
were dried using the HAD technique. The 
experiments were carried out in a pilot plant 
hot-air drier (UOP 8 Tray dryer, Armfield, 
UK, equipped with automatic data recording 
system and temperature and airflow velocity 
controller units) at 45, 55, and 65+41°C and the 
airflow rate of 1 m/s (Doymaz 2012; 
Lemus-Mondaca et al., 2015). Moisture loss 
was calculated by measuring the mass loss of 
the samples at 15 min intervals (based on 
preliminary experiments) by a precision 
balance with an accuracy of +0.01 g. Moisture 
content data were recorded throughout the 
drying experiments using a data logger 
connected to a PC. The experiments were 
continued until reaching a final moisture 
content of 10% (wet basis). 


Modeling of drying curves 


Eighteen different empirical and semi- 
empirical models were used to evaluate the 
kinetics of moisture loss during HAD of lemon 
verbena leaves (Table 1) (Ertekin and Heybeli 
2014). The model parameters or drying 
constants (a, b, c, g, h, k, and n) were 
estimated by applying non-linear regression 
analysis using MATLAB software (Version 
8.1.0.604 R2013a, The Math works, Inc., 
USA). The coefficient of determination (R’), 
adjusted R’, root mean squared error (RMSE), 
and sum of squared error (SSE) were used to 
evaluate the goodness of fit in order to select 
the suitable model(s) to predict the drying 
kinetics. These statistical criteria are as 
follows (Lemus-Mondaca et al., 2015): 


N 
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1 2 
RMSE = (ze (MR a, ;- MR ~?| (2) 


1 N 
SSE = Re (MR oxi MR pre)” (3) 
i=l 


Where MRexp,i is the i experimental 


moisture ratio, MRpre,i Shows the ith predicted 
moisture ratio, MRexp stands for the average 


experimental moisture ratio, and N denotes the 
number of observations or the number of data 
values. 


Table 1. Kinetic models used to describe the drying of lemon verbena leaves" 


Model number Model equation 


Model name 


1 MR = expC kt) 


Lewis (Newton) 


2 MR =exp(—kt") Page 
3 MR =exp(—(kt)") Modified Page-I 
4 MR = aexp(—kt) Henderson & Pabis 
5 MR = aexp(—kt") Modified Page-II 
6 MR =aexp(—kt") + bt Midilli et al. 
7 MR = aexp(—kt)" +b Demir et al. 
8 MR =(a—b)exp(—kt") Weibull distribution-I 
9 MR = (a— b) exp(—(kt)") Weibull distribution-II 
10 MR= exp) Aghlasho 
11 ———— Logistic 
1+ bexp(kt) 
12 MR = aexp(—kt) + bexp(— gt) Two-term 
13 MR = aexp(— kt") + bexp(- gt”) Hii et al. 
14 MR = aexp(— kt) + (1 — a) exp(— kat) Two-term exponential 
15 MR = aexp(—kt) +(1—a)exp(-bt) Modified two-term exponential 
16 MR = aexp(— kt) + (1 — a) exp(— kbt) Diffusion approximation 
17 MR = (1—a)exp(—kt) + (1— b) exp(—gt") +c Naghavi et al. (present study) 
18 MR = aexp(—kt) + bexp(— gt) + cexp(—ht) Modified Henderson & Pabis 


*All models (except model-17) are available in the paper published by Ertekin and Heybeli (2014). 
a, b, c, g, h, k, and n are model parameters (empirical constants). 


Determination of the _ effective moisture 
diffusivity 

The Fick’s law-based model (Eq. 4) is often 
used to determine the effective moisture 


diffusivity (Derr) of different food materials. 
2 
eM =) ett o™ (4) 
Ot Ox 

The following initial and boundary 
conditions can be considered (Doymaz 2012): 

t=0, O<x<L, M=Mo 

t>0O, x=L, M=Me 


t>0O, x=0, ete 


In the present study, the analytical solution 
of Fick’s second law for an infinite slab (Eq. 
5) was applied to calculate Der (Erbay and 
Icier 2010; Doymaz 2012): 


mre-MM_ 8S 1 { onset] 
(5) 


M-M, 2 Qn+1y 4V 
where M is the moisture content (dry basis), 


Detf represents the effective moisture 
diffusivity (m7/s), L is the half thickness of the 
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slab (m), t stands for the time (s), MR is the 
moisture ratio (dimensionless), M: shows the 
moisture content at any time (kg water/kg dry 
solid), Mo is the initial moisture content (kg 
water/kg dry solid), Me denotes the 
equilibrium moisture content (kg water/kg dry 
solid), and n is the number of the terms taken 
into consideration. 

This method has also been used previously 
by several researchers (Doymaz 2012; Aral 
and Bese 2016) and is based on_ the 
assumptions that shrinkage is negligible, Derr 
remains constant and moisture loss occurs 
through the diffusion phenomenon (Crank 
1975). For long drying times (M-=0), the use 
of one-term approximation (n=1) to the series 
summation is reasonable and Eq. 5 reduces to 
(Doymaz 2012): 


M, ee _@Dagt (6) 
M, « Aap 


By taking the natural logarithm, Eq. 6 can 
be further simplified to (Doymaz 2012; Oberoi 
and Sogi 2015): 


In(MR)= in( 8] -(* Pet) (7) 
1 


A plot of the experimental data in terms of 
In (MR) versus time gives a straight line with a 
slope of ko (Oberoi and Sogi 2015): 
-1’D, 
ky > 4I2 = (8) 


MR = 


Determination of the activation energy 


The Arrhenius model describes the 
relationship between Derr, drying temperature 
(T), and the activation energy (Ea). Therefore, 
for quantifying Ea and investigating the effect 
of temperature on Derr, the Arrhenius type 
equation (Eq. 9) was employed (Doymaz 
2012): 


E 
Deg = b,ev{- RT (9) 


Where E, denotes the activation energy 
(kJ/mol), R is the universal gas constant [8.314 
J/(mol K)], T is the absolute air temperature 
(K), and Do is the pre-exponential factor 
(constant) (m?/s). This approach has been used 


by a number of researchers (Doymaz 2012; 
Aral and Bese 2016). 


Determination of the rehydration ratio 


Measurement of the rehydration ratio (RR) 
of dried leaves was carried out according to 
the method described by Doymaz et al. (2015) 
and Nozad et al. (2016). Based on this method, 
5 g of the dehydrated samples were poured 
into a glass beaker (750 mL) containing 500 
mL of distilled water (25°C) and kept for 24 h. 
Next, the leaves were removed from the 
beaker and their surface water was blotted up 
using a tissue paper. Finally, the weight of the 
resulted sample was measured precisely using 
a digital balance. In all cases, the tests were 
triplicated for each sample and the mean 
values of the three replications + standard 
deviation were reported. The RR calculation 
was carried out using Eq. 10 as employed by 
Doymaz et al. (2015), Nozad et al. (2016), and 
Salarikia et al. (2016): 

= oe W, 
RR = W, (10) 

Where RR denotes the rehydration ratio [kg 
water/kg dry matter (DM)], W1 is the weight 
of the dried leaves (kg), and W2 represents the 
weight of the rehydrated leaves (kg). 


Color measurement 


The color changes of leaves (fresh and hot- 
air dried) were quantified using image 
processing in MATLAB (Version 8.1.0.604 
R2013a, The Math works, Inc., USA) (Nozad 
et al. 2016). The color test instrument was 
designed and constructed in the Department of 
Agricultural Machinery Engineering, 
University of Tabriz, Tabriz, Iran. It consists 
of a chamber with a trapezoidal cross section 
that was equipped by two Des (daylight) lamps 
as the light source for illumination of sample. 
At first, a sample was put in the chamber. 
After zooming the lens and focusing, the 
images were taken by camera. A digital 
camera (Nikon, D3200, Japan) was used to 
capture images from leaf surfaces. The camera 
calibration was performed prior to each drying 
experiment. 
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In each experimental run, the color of the 
leaves (10 fresh and 10 dried samples) was 
measured as L*, a* and b* values, which 
known as Hunter parameters. It is well known 
that the L* value represents the degree of 
lightness/darkness, a* stands for the degree of 
redness/greenness, and b* shows the degree of 
yellowness/blueness. Changes in the color of 
the leaves were calculated as follows (Nozad 
et al., 2016; Salarikia et al., 2016): 

AB=VAL* +Aa® +Ab® (11) 

Where AE denotes the total color change of 
leaves and AL*, Aa*, and Ab* represents the 
difference between the color parameters of 
initial samples (Lo*, ao*, and bo*) and final 
dried leaves (L*, a*, and b*). 


Determination of the essential oil content 


The essential oil was extracted from the 
leaves using a flask connected to Clevenger 
hydro-distillation apparatus (Nozad eft al., 
2016). Based on this method, a given amount 
of the samples (30 g) were put into a round- 
bottomed distillation flask filled with a given 
amount of distilled water (250 mL). Then, the 

1 


heating was performed for 3 h and the distilled 
essential oil collected in the side arm was 
separated. The data of essential oil (%) were 
expressed on the basis of dry matter weight. 


Statistical analysis 


The experimental data were analyzed 
statistically by the analysis of variance 
(ANOVA) using Minitab statistical software 
(Minitab Release 14, Minitab Inc., USA). The 
significant difference between the means was 
determined using Tukey’s honestly significant 
difference (HSD) test at the significance level 
of 5% (p<0.05). The data were expressed as 
the mean + standard deviation and all 
experiments were carried out in triplicate. 


Results and discussion 
Drying kinetics 


The effects of the drying temperature (45, 
55, and 65°C) on the dimensionless moisture 
ratios of lemon verbena leaves are illustrated 
in Fig. 1. 
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Fig. 1. Moisture ratios (dimensionless) as a function of the drying time at different temperatures during hot-air drying of 
lemon verbena leaves 


It can be seen that the moisture ratio 
declined quickly in the initial period of HAD 
(almost up to 75-90 min) and was 
subsequently followed by a gradual non-linear 
decrease (almost exponential) with an increase 
in the process time. This result was similar to 
the findings of other researchers (Doymaz 


2012; Lemus-Mondaca et al., 2015; Said et al., 
2015; Aral and Bese 2016). 

Fig. 1 also shows that all drying 
temperatures exhibited a relatively similar 
behavior for the dried samples. Moreover, by 
increasing the drying temperature, the 
moisture ratios and consequently the drying 
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kinetics were altered. Based on ANOVA 
results, the drying time was_ significantly 
(p<0.05) lower in the case of the samples dried 
at 65°C than those dehydrated at lower 
temperatures (45 or 55°C). The total drying 
time was 420, 345, and 285 min for the leaves 
dried at 45, 55, and 65°C, respectively. This 
result indicates that a temperature increase of 
20°C (1.e. from 45 to 65°C) caused a reduction 
of approximately 135 min in the total drying 
time (p<0.05). This might be due to the 
increased vapor pressure in lemon verbena 
leaves at higher drying temperatures, which in 
turn results in a faster moisture loss from the 
samples and thus, a shorter drying time (Aral 
and Bese 2016). Other studies have reported 
similar findings (Ertekin and Heybeli 2014; 
Lemus-Mondaca et al., 2015; Said et al., 2015; 
Aral and Bese 2016). 


Modeling of the drying curves 


In the literature, several empirical and semi- 
empirical models have been applied to 
describe the drying kinetics of food materials 
and plants under different drying conditions, 
which are based on the dimensionless moisture 
ratio as a function of drying time. In the 
present work, the experimental data of 
moisture loss in lemon verbena leaves during 
HAD at different drying temperatures (45, 55, 
and 65°C) were fitted to eighteen models 
summarized in Table 1 (Ertekin and Heybeli 
2014), which allows of using and comparing 
different model correlations based on the 
fitting criteria (R’, RMSE, and SSE) on a 
defined set of experimental data 


Table 2. Statistical analysis of different kinetic models for different drying temperatures 


Model T 2 Adjusted Model T 2 Adjusted 
aarabee (°C) R R? RMSE SSE aiAbed (°C) R R? RMSE SSE 

1 45 0.9875 0.9875 0.0283 0.0225 10 45 0.9977 0.9976 0.0124 0.0041 
55 0.9870 0.9870 0.0294 0.0199 55 0.9981 0.9980 0.0116 0.0029 
65 0.9907 0.9907 0.0256 0.0125 65 0.9981 0.9980 0.0120 0.0026 

2 45 0.9995 0.9995 0.0056 0.0008 11 45 0.9744 0.9735, 0.0413 0.04611 
55 0.9996 0.9996 0.0050 0.0005 55 0.9723 0.9710 0.0439 0.04233 
65 0.9997 0.9997 0.0049 0.0004 65 0.9766 0.9753 0.0417 0.03132 

3 45 0.9995 0.9995 0.0056 0.0008 12 45 0.9998 0.9998 0.0034 0.00028 
35 0.9996 0.9996 0.0050 0.0005 55 0.9998 0.9999 0.0032 0.00019 
65 0.9997 0.9997 0.0049 0.0004 65 0.9999 0.9998 0.0033 0.00017 

4 45 0.9927 0.9924 0.0221 0.0131 13 45 0.9999 0.9999 0.0030 0.00021 
55 0.9916 0.9912 0.0241 0.0128 55 0.9997 0.9996 0.0051 0.00049 
65 0.9934 0.9931 0.0221 0.0088 65 0.9997 0.9996 0.0053 0.00042 

| 45 0.9995 0.9995 0.0057 0.0008 14 45 0.9996 0.9996 0.0048 0.00063 
55 0.9996 0.9996 0.0051 0.0005 55 0.9995 0.9994 0.0061 0.00083 
65 0.9997 0.9996 0.0050 0.0004 65 0.9998 0.9998 0.0041 0.00030 

6 45 0.9997 0.9997 0.0044 0.0005 15 45 0.9998 0.9998 0.0033 0.00029 
55 0.9997 0.9997 0.0046 0.0004 55 0.9999 0.9998 0.0032 0.00027 
65 0.9998 0.9997 0.0044 0.0003 65 0.9999 0.9998 0.0031 0.00026 

7 45 0.9997 0.9997 0.0047 0.0005 16 45 0.9998 0.9998 0.0033 0.00030 
55 0.9997 0.9996 0.0048 0.0005 55 0.9999 0.9998 0.0030 0.00025 
65 0.9998 0.9997 0.0044 0.0003 65 0.9999 0.9999 0.0032 0.00017 

8 45 0.9995 0.9995 0.0058 0.0008 17 45 1 1 0.0007 0.000011 
35 0.9996 0.9996 0.0052 0.0005 55 1 0.9999 0.0013 0.000031 
65 0.9997 0.9996 0.0052 0.0004 65 0.9999 0.9999 0.0020 0.000056 

9 45 0.9995 0.9995 0.0058 0.0008 18 45 0.9998 0.9998 0.0035 0.00028 
55 0.9996 0.9996 0.0052 0.0005 55 0.9999 0.9998 0.0035 0.00022 
65 0.9997 0.9996 0.0052 0.0004 65 0.9997 0.9997 0.0049 0.00034 


Table 3. Estimated coefficients of different kinetics models at different drying temperatures’ 
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Model number _ T (°C) Coefficients Model number — T (°C) Coefficients 
1 45 a=0.0102 10 45. a=0.0124 b=0.0014 
55 a =0.0130 55. a=0.0159 b=0.0019 
65 a =0.0177 65 a=0.0212 b=0.0022 
2 45 a=0.0230 n=0.8315 u 45 a=1.7520 k=0.0128 
55 a=0.0294 n=0.8213 55 a=17650 k=0.0164 
65 a=0.0360  n=0.8346 65 a=18030 k=0.0229 
a =0.0107 a=0.7404 g = 0.0350 
2 - n=0.8315 ” 45h = 0.2604 k= 0.0078 
= a =0.0137 ss 2=0.7025  g = 0.0402 
n=0.8213 b=0.2977 k=0.0096 
- a = 0.0186 65 80.7515 g = 0.0625 
n= 0.8346 b=0.2490 k=0.0138 
a = 0.6834 
= 0.9339 k = 0.0051 
4 45 a 1B 45. b=0.3152 
k = 0.0095 gzo.o2g3 2 1.0690 
a= 1.0290 
a =0.9370 k = 0.0338 
55 ; i hep 0090 
k =0.0121 ai 220.7970 
a = 0.9833 
a = 0.9500 k = 0.0334 
i) k = 0.0168 a : sey 20.8485 
a= 1.0020 : a = 0.2325 
5 45 ae ge = 0.8298 14 doh Soe 
a= 1.0020 a = 0.2449 
22 coo. “BEOdES 55 k= 0.0417 
a = 0.9993 7 a =0.2307 
os) cepessg RRS 65k =0.0612 
a=1.0060  k=0.0259 a = 0.2599 
6 45h =-3.089x105 n=0.8041 b > Gpeoaare” SOO 
a=1.0040 k =0.0317 a =0.2977 
55h =-2.328x105 n= 0.8036 S: penjoe <=! 
a=1.0010  k =0.0380 a=02487 || 
65h =-2.562x10° n=0.8199 OF Reise eet 
a=1.0210  k=0.0104 a=02598 | _ 
‘ #2 b=-0.0150 n=0.8015 i oP egsag. 6 ee 
a=1012 k=0.0135 a=02977 | | 
3h =-0.0084 nn = 0.8030 2 _ yegoaug = OHO 
a=1.0100  k=0.0183 a =0.2490 
© b=-0.0087 n=0.8154 O° - pegooia: “SS O0ee3 
= 7 a=0.6936 g=0.0114 
8 45 a ae pee 7 45. b=0.2731 k=0.0256 
=a), = ¢ = 0.0333 n=0.9001 
a=1.1360  k=0.0298 Rae rece ae yr 
55 Soy cera 55. b=1.9780  k=0.0153 
=i =: 6=06156 n=-0.4296 
: 5 a=03155 g=0.0001 
65 Poa! ones 65 b=0.6998 k =0.0352 
aie ae c=0.0144 n=1.8710 
: : a=0.7466 g=0.0373 
9 eS SS Ben ae es 18 45. b=0.2643 h=10.300 
=-7. =) ¢=-0.0109 k =0.0079 
7 : a=-0.1141 g=0.0495 
oe 55. b=0.4036 hh =0.0096 
= is c=0.7106 k=0.0802 
a=0.5256 k=0.0186 Be Us eee 
Gi Fae aa 65 b=-1424  h=0.2176 
= ae c=16190 k=0.0144 


* a, b,c, g, h, k, and n are model parameters (empirical constants). 
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Fig. 2. Comparison of the experimental (exp.) and predicted (pred.) moisture ratios for lemon verbena leaves dehydrated at 
different drying temperatures 


Table 2 indicates that for the most of tested 
models, R? and adjusted R* values were 
higher than 0.999, and RMSE and SSE values 
were between 0.0007-0.0413 and 0.00001 1- 
0.04611, respectively. 

Genrally, the closer the experimental and 
predicted moisture ratios, the better they 
explain the adequacy of the regression model. 
As expected, the models with a larger number 
of coefficients (models 12, 13, 15-18 in Table 
2) had higher R? and lower RMSE and SSE 
values. The data in Table 2 suggested that the 
employed models were suitable to describe 
the drying behavior of lemon verbena leaves. 
However, for plotting the predicted moisture 
ratios against the drying time, only the model 
presented in this research (model-17 in Table 
1) was fitted to the experimental data (Fig. 2) 
which had greater R* (0.9999-1) and adjusted 
R? (0.9999-1) values and lower RMSE 
(0.0007-0.0020) and SSE  (0.000011- 
0.000056) values than the other 17 models 
(Table 2). The models coefficients (constants) 
obtained at different drying temperatures are 
represented in Table 3. By increasing drying 
temperature, different coefficients did not 
follow a similar trend. As can be obtained 
from Table 3, drying rate constant (k) in the 
studied drying models increased with 
increasing drying temperatures. Thus, it may 
be assumed that this kinetic parameter would 


be directly proportional to drying temperature 
(Lemus-Mondaca et al., 2015). Similar results 
were reported by other investigators (Vega- 
Galvez et al., 2012; Lemus-Mondaca et al., 
2015). Furthermore, there was no clear trend 
on the effect of drying temperature on the 
other constants (a, b, c, g, h, and n). 


Rehydration ratio (RR) 

Comparison of the RR results for the 
samples dried at different drying temperatures 
are presented in Fig. 3. 

It can be seen that there was no significant 

change (p>0.05) in the RR of the dried leaves 
when the drying temperature changed from 45 
to 55°C. However, RR reduced significantly 
(p<0.05) with an increase in the air 
temperature from 45 to 65°C. This was 
attributed to the fact that HAD at higher 
temperatures resulted in a change in the 
structural features of the leaves (such as tissue 
collapse, development of a surface hard layer 
and volumetric shrinkage), which in turn can 
cause significant damages to the textural 
quality of the dried samples and therefore a 
decrease in RR (Doymaz 2012; Nozad et al., 
2016; Salarikia et al., 2016). The highest RR 
value of the dried samples (78.24%) was 
observed for the leaves dried at 45°C, 
followed by those dehydrated at 55 (76.12%) 
and 65°C (73.36%). These results are in 
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agreement with the results reported by f 
Doymaz (2012), who stated that the higher 
the drying temperature (40, 50, and 60 °C) 
resulted in a lower RR of the grape leaves, 
concluding that higher drying temperatures 
led to greater changes in the structural 
attributes and thus, lower RR values. Nozad et 
al., (2016) also observed that in HAD of 

0.8 


0.78 
0.76 
0.74 


0.72 


Rehydration ratio (kg water/ kg DM) 


0.7 


45 


spearmint (Mentha spicata L.) leaves, an 
increase in the air temperature from 30 to 
50°C had a considerable decreasing effect on 
RR. Similar findings were reported by other 
researchers (Jangam ef al., 2008; Vega- 

Galvez et al., 2012; Salarikia et al., 2016). 


ye) 65 


Drying temperature (°C) 


Fig. 3. Rehydration ratio of lemon verbena leaves at different drying temperatures. Error bars show one standard deviation 
from the mean and means with different letters are significantly different (p<0.05). 


Effective moisture diffusivity (Desr) 


In the present study, the value of Derr at 45, 
55, and 65°C was estimated by plotting 
In(MR) versus the drying time (Eq. 4) (Oberoi 
and Sogi 2015). The slope of _ the 
corresponding line represents the value of Derr. 
The Derr values varied from 1.140 x 10°'° to 
2.280 x 10° m?/s (Table 4). Other researchers 
found similar values for the Derr of dried foods 
(in particular leaf materials), which were in 
general in the range of 10°” to 10° m/’/s 
(Zogzas et al., 1996). Derr values have been 
determined for other herbs as follows: mint 
leaves 3.067 x 10° to 1.941 x 10° m’/s 
(Doymaz 2006), spinach leaves 6.590 x 10°! 


to 1.927 x 10° m/s (Doymaz 2009), nettle 
leaves 1.744 x 10° to 4.992 x 10° m*/s (Kaya 
and Aydin 2009), mint leaves 0.965 x 10°'! to 
1.190 x 10°!! m/s (Therdthai and Zhou 2009), 
olive leaves 1.054 x 10° to 4.973 x 10° m?/s 
(Erbay and Icier 2010), grape leaves 4.13 x 10° 
10 to 1.83 x 10° m/s (Doymaz 2012), kaffir 
lime leaves 2.61 x 10°!! to 9.24 x 10°!! m?/s 
(Tasirin et al., 2014), stevia leaves 4.67 x 10°!° 
to 14.90 x 10°!° m?/s (Lemus-Mondaca et al., 
2015), and wild edible plant (Allium roseum) 
leaves 2.55 x 10°! to 8.83 x 10°!* m’/s (Said et 
al., 2015). 


Table 4. The values of effective moisture diffusivity (Dest) for lemon verbena leaves at different drying temperatures 


Drying temperature (°C) 


Deer (m/s) 


45 
55 
65 


1.140 x 10°! ¢ 
1.710 x 1071 
2.280 x 10° 


Different letters in the same column indicate significant differences (p<0.05) 


As can be realized from Table 4, significant 


differences (p<0.05) of the Derr values were 
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observed between the leaves dried at different 
drying temperatures. This indicates that hot-air 
temperature has a considerable effect on Dest 
during HAD of plants and foods, as reported 
by numerous authors (Erbay and Icier 2010; 
Doymaz 2012; Lemus-Mondaca ef al., 2015; 
Said et al., 2015; Aral and Bese 2016). This 
may be related to the higher thermal energy 
transferring to the leaves at higher drying 
temperatures, which subsequently results in an 
increase in the kinetic energy of the water 
molecules (Aral and Bese 2016). A higher Der 
value indicates the increasing rate of moisture 
loss with the rise of drying temperature (Fig. 
1). 


Activation energy (Ea) 


The value of E, in the drying of lemon 
verbena leaves was estimated from the slope 
of the linearized Arrhenius equation (Eq. 9) 
(Doymaz 2012) and was found to be 31.04 
kJ/mol. The Ea value obtained in this research 
was in the range that reported for other 
aromatic plants and fruits. Experimentally- 
determined E, values have been reported by 
several researchers, for example Ahmed et al. 
(2001) for coriander leaves (26.50 kJ/mol in 
the temperature range of 45-65°C), Doymaz 
(2006) for mint leaves (62.96 kJ/mol in the 
temperature range of 35-60°C), Doymaz 
(2009) for spinach leaves (34.35 kJ/mol in the 
temperature range of 50-80°C), Kaya and 
Aydin (2009) for nettle leaves (79.873- 
109.003 kJ/mol in the temperature range of 35- 
55°C and at airflow rates of 0.2-0.6 m/s), 
Erbay and Icier (2010) for olive leaves (60.97 
kJ/mol in the temperature range of 50-70°C), 
Doymaz (2012) for grape leaves (64.56 kJ/mol 
in the temperature range of 40-60°C), Lemus- 
Mondaca et al. (2015) for stevia leaves (38.78 
kJ/mol in the temperature range of 30-80°C), 
and Said et al. (2015) for wild edible plant 
(Allium roseum) leaves (46.80-52.68 kJ/mol in 
the temperature range of 30-80°C and 1 and 1.5 
m/s airflow velocity). 

It has been reported that the value of Ea is 
influenced by several factors, including the 
drying air temperature, the moisture content of 
food or herb, and variations in the Derr value 


with the drying temperature (Aghbashlo et al., 
2008), which makes it difficult to compare the 
E, values for different foods and _ herbs 
dehydrated at different process conditions. 
However, from numerous conducted studies 
on the calculation of the E, value, it can be 
concluded that long dehydration time, high 
initial moisture content, remarkable variation 
in the Derr value with the drying temperature 
(at constant airflow rate) (Aghbashlo et al., 
2008) or with both temperature and airflow 
rate (Erbay and Icier 2010), low hot-air flow 
rate, low drying temperature, and textural 
changes in the sample due to the percentage of 
shrinkage and tissue collapse, are all the 
reasons for a considerable increase in the E, 
value. 


Color measurement 

The value of color change (AE) represents 
the degree of total color change in dehydrated 
leaves compared to the color of fresh samples. 
The lower AE the better the quality of dried 
leaves (Salarikia et al., 2016). The color of 
aromatic plants and herbs are very sensitive to 
heat damage during HAD. Fig. 4 shows the 
value of AE for dried samples. It can be seen 
that the value of AE increase significantly 
(p<0.05) with increasing of hot-air temperature 
from 45 to 65°C. Chlorophyll a and 
chlorophyll b are responsible for natural green 
color of leaves. The increase in AE is due to 
the increase in substitution of magnesium with 
hydrogen in chlorophyll with drying 
temperature. Under this condition, 
chlorophylls are converted to pheophytins 
(Therdthai and Zhou 2009). This finding 
confirmed the previous observations obtained 
by Therdthai and Zhou (2009) for mint leaves 


(Mentha cordifolia Opiz ex _ Fresen), 
Chenarbon et al. (2012) for St. John's 
wort (Hypericum —_perforatum _L.) leaves, 


Akbudak and Akbudak (2013) for parsley, and 
Salarikia et al. (2016) for peppermint leaves. 


Determination of the essential oil content 
Essential oil content of lemon verbena 

leaves before (fresh sample) and after drying is 

shown in Fig. 5. The fresh samples had the 
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highest value of essential oil (1.10%) between in an increase in rupture of oil glands and as a 
the treatments. Furthermore, no significant consequence, rapid evaporation or higher loss 
difference (p>0.05) of essential oil content was of volatile components (Argyropoulos and 
observed between the leaves dried at different Miiller 2014). This result is consistent with 
temperatures, while significant difference those obtained in previous studies, for example 
(p<0.05) of essential oil content was observed Laurus nobilis L. leaves (Sellami et al., 2011), 
between fresh leaves and samples dehydrated Thymys daenensis subsp. daenensis. Celak 
at 55 and 65°C. The reduction in essential oil leaves (Rahimmalek and Goli, 2013), Lemon 
content with increasing drying temperature verbena (Shahhoseini et al., 2013), and 


(p>0.05) might be explained by the fact that peppermint leaves (Salarikia et al., 2016). 
the relatively high temperature of hot-air result 


16 4 a 7 


Color change (AE) 


Drying temperature (°C) 


Fig. 4. Color change in lemon verbena leaves at different drying temperatures. Error bars show one standard deviation from 
the mean and means with different letters are significantly different (p<0.05). 
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Fig. 5. Essential oil content of lemon verbena leaves at different drying temperatures .Error bars show one standard 
deviation from the mean and means with different letters are significantly different (p<0.05). 


Conclusions some selected properties of lemon verbena 


This study was focused on investigation of leaves (moisture loss, Deft, color change, 
essential oil content, and RR) during HAD at 
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different drying temperatures as well as the 
empirical modeling of the drying kinetics. A 
shorter drying time and a higher Derr value 
were observed (p<0.05) with an increase in the 
drying temperature from 45 to 65°C. The 
values of Derr ranged from 1.140 x 10°'° to 
2.280 x 10° m?/s and the Ea value was found 
to be 31.04 kJ/mol, all of which were in 
agreement with the results reported by other 
investigators in the literature. Our results also 
showed that the percentage of RR _ was 
significantly (p<0.05) affected by the air 
temperature and its maximum (78.24%) and 
minimum (73.36%) values were attained at 
45°C and 65°C, respectively. Essential oil 


temperatures was not significant (p>0.05) with 
respect to each other but was significantly 
(p<0.05) different from the fresh samples. 
Furthermore, the value of AE increased 
significantly (p<0.05) with increasing of hot- 
air temperature from 45 to 65°C and ranged 
from 12.24 to 14.92, respectively. The results 
of modeling (R?>0.99 and low RMSE and SSE 
values for most of the tested models) indicated 
a good fit to the experimental data of moisture 
ratio. Among these, the model proposed in the 
present study had a better goodness of fit (with 
an adjusted R* >0.999 and the lowest RMSE 
and SSE) and was considered as the best 
model. 


content of the samples dried at different drying 
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